Microscopic theory of temperature-dependent magnetoelectric effect in Cr 2 0,3 
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We calculate the temperature-dependent magnetoelectric response of Cr2 03 from first principles. 
The form of the dominant magnetoelectric coupling is determined using symmetry arguments, its 
strength is found using ab initio methods, and the temperature dependence of the response is 
obtained from Monte Carlo simulations. The quantitative agreement of our results with experiment 
shows that the strong temperature dependence of the magnetoelectric effect in Cr2C>3 results from 
non-relativistic exchange interactions and spin fluctuations. 
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Introduction: Recent progress in understanding the sub- 
class of multifcrroics in which improper ferroelectric po- 
larizations are induced by non-centrosymmetric magnetic 
orderings has led to a clarification of the microscopic ori- 
gins for magnetoelectric coupling. 1 - In particular two dis- 
tinct coupling mechanisms have been identified. The first 
arises from relativistic effects linking electron spin and or- 
bital momentum, resulting in the antisymmetric Si x S2 
interaction between spins of magnetic ions. The depen- 
dence of the strength of this Dzyaloshinksii-Moriya in- 
teraction on polar displacements of ions makes magnets 
with non-collinear spiral orders ferroelectric In the 
second mechanism, polar deformations of the lattice are 
induced by Heisenberg spin exchange interactions Si • S2, 
originating from the Fermi statistics of electrons*^ This 
non-relativistic mechanism can give rise to stronger mag- 
netoelectric couplings than those resulting from relativis- 
tic effects, which tend to be relatively weak in 3d tran- 
sition metal compounds, and is not restricted to non- 
collinear spin arrangements. Indeed, the electric polar- 
izations of the orthorombic manganite Yi-^Lu^MnOa 8 
and orthoferrite GdFe03,? which have collinear antifer- 
romagnetic spin orderings and polarization arising from 
the Heisenberg mechanism, exceed the largest polariza- 
tions observed in spiral multiferroics by one order of mag- 
nitude. 

In this work we show that, in addition to causing multi- 
ferroic behavior, the relativistic and Heisenberg exchange 
mechanisms can both give rise to the linear magneto- 
electric effect, in which an applied magnetic field in- 
duces an electric polarization proportional to the field 
and, conversely, magnetization is induced by an applied 
electric field (for a review see Reffioh. Furthermore, 
the two mechanisms can co-exist in the same material 
and dominate in different temperature regimes. Using a 
combination of first-principles density functional theory 
and Monte Carlo methods we calculate the temperature- 
dependent magnetoelectric response of the prototype 
magnetoelectric material, chromium sesquioxide, Cr203. 
We show that the strong finite-temperature magnetoelec- 
tric response originates from the Heisenberg exchange 
mechanism combined with thermal spin fluctuations, 
whereas at zero temperature, where spin fluctuations 




FIG. 1: (Color online) Rhombohedral unit cell of Cr203 with 
the unit vectors a,, i = 1,2,3, containing four magnetic Cr 
ions. Panel (a) shows the actual antiferromagnetic spin or- 
dering in Cr2 03, while panel (b) shows the spin ordering im- 
posed in our first-principle calculations to induce an electric 
polarization along the trigonal z axis. 



vanish, the observed weak response arises from relativis- 
tic effects. Significantly, Cr203 is a collinear antifcrro- 
magnet, and so our work extends the recent suggestion 
that Heisenberg exchange interactions can give rise to 
a strong magnetoelectric effect in frustrated magnets in 
which spins are forced to form non-collinear orders with 
nonzero toroidal or monopole momentsj 11 ' 12 
Magnetoelectric coupling: The magnetoelectric ef- 
fect in Ci - 203 was predicted phenomenologically by 
Dzyaloshinksii 1 ? and measured by Astro\*i^ shortly after 
the theoretical prediction. Cr203 has four Cr 3+ ions with 
spin S — 3/2 in the rhombohedral unit cell. Below T/y 
= 307K it shows collinear tiTi s P m ordering along the 
trigonal z axis 15 [see Fig. [T](a)]. This antiferromagnetic 
ordering preserves the three-fold and two-fold symmetry 
axes of the paramagnetic phase, but it breaks inversion 
symmetry /, reducing it to I' (inversion combined with 
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time reversal), which allows for two independent magne- 
toelectric coupling terms in the free energy^: 

Fme = -9\\G Z E Z H Z - g ± G z {E X H X + E y H y ) , (1) 

where G z — (S* — <Sf + 5| — S£) is the antiferromag- 
netic order parameter. The two magnetoelectric coeffi- 
cients, ay = g\\G z and aj_ = g±_G z , show very different 
temperature dependences; 14 While a± varies only in the 
immediate vicinity of the transition temperature, T/v, the 
coefficient au is very strongly temperature dependent. 
a|| reaches its maximum at T max ~ 260K, below which 
it steeply decreases and changes sign at ~ 100K. While 
at T max the magnitude of ck | is one order of magnitude 
larger than |aj_|, at low temperatures a|| < |ot_i_ j . 

Early measurements of the temperature dependence 
of the magnetoelectric coefficients* 1 ^— as well as re- 
cent first-principles calculations, 1 - show that the rela- 
tively weak magnetoelectric effects at low temperature 
result from relativistic interactions. In what follows, we 
demonstrate that the much stronger response at elevated 
temperatures originates from Heisenberg exchange. 

Phenomenologically, the coupling of the electric polar- 
ization P z along the trigonal axis to spins, resulting from 
the dependence of spin-exchange constants on polar lat- 
tice distortions, is described by 

P 2 = A(S 1 -S 3 -S 2 -S 4 ), (2) 

where S Q with a = 1, 2, 3, 4 denotes the sublattice mag- 
netization and A is the coupling strength that we will de- 
termine from ab initio calculations. Here we take into ac- 
count the fact that the exchange-driven polarization can 
only depend on scalar products of the magnetizations. 
The combination of scalar products in the right-hand side 
of Eq. ([2]) transforms as P z , which changes sign under C2 
and I. The combination of scalar products that lead to 
this property can be seen by inspection of Table HI which 
shows how the four inequivalent magnetic sites transform 
under the symmetry operations of Cr 2 3 . 
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TABLE I: Transformation of four independent Cr sites with 
the fractional coordinates ri = (u, it, it), r 2 = (1/2 — u, 1/2 — 
U, 1/2 - u), r 3 = (1/2 + u, 1/2 + u, 1/2 + it), and r 4 = 
(1 — it, 1 — it, 1 — it), where it ~ 0.153, under the genera- 
tors of space group R3c: the 120°-rotation around the z axis, 
C3 = (x3,Xi,X2), the 180°-rotation around the axis orthogo- 
nal to the ^-direction, C2 = (1/2 — x%, 1/2 — xi, 1/2 — £3) and 
inversion I = (1 — xi, 1 — x 2 , 1 — X3). Here, c = ai + a 2 + sls, 
where ai (i = 1, 2, 3) are the rhombohedral unit vectors. 

Equation © is clearly appropriate for describing elec- 
tric polarization induced by spin ordering in a multi- 
ferroic material. In addition, it applies to the linear 



magnetoelectric effect, as can be seen from the follow- 
ing heuristic argument: In an applied magnetic field H z 
the average value of spin on the sublattice a changes 
by (SS^) oc X\\H Z , where x\\ is the longitudinal mag- 
netic susceptibility. Equation ([2]) then gives P z oc 
Ax|| {SI -S$+SI- SI) H z oc G Z H Z , consistent with 
pz = _^kp = g l{ G z H z obtained from Eq.©. 

Equation (2) is only meaningful within the mean field 
approach. To account for effects of spin fluctuations on 
the magnetoelectric response of Cr2C>3 we will use the 
microscopic expression for the exchange- driven polariza- 
tion in terms of scalar products of Cr spins (rather than 
the sublattice magnetizations), which has the form 

A 6 

pZ = 6iV X! ( Sl >3 ' S 3,j-b™ ~ S 2,J ' S 4J-b„ ) • (3) 
j n=l 

Here j labels unit cells, TV is the total number of unit 
cells, and b : = a l7 b 2 = a 2 , b 3 = a 3 , b 4 = a x + a 2 , 
= a 2 + a 3 b 6 = a 3 + a 4 (a, being the rhombohe- 
dral unit vectors). Remarkably, the interaction between 
the fourth nearest-neighbor Cr ions separated by the dis- 
tance ri ! ( 3 _ Q1 ) = r li(3 _ ai _ a2 ) = 3.65A turns out to give 
rise to the magnetoelectric effect in Cr 2 3 . The shorter- 
range exchange interactions do not couple the sublattices 
1 and 3 (or 2 and 4) and, therefore, do not contribute 
to P z , while other interactions between these sublattices 
correspond to much longer exchange paths and are neg- 
ligibly small. This allows us to accurately pinpoint the 
microscopic origin of the strong magnetoelectric effect in 
Cr 2 3 . 

Importantly, Eqs.© and §S§ apply to any four- 
sublattice spin ordering in Cr 2 3 . The t4-T4- s P m or ~ 
dering realized in the low-temperature ground-state of 
Cr 2 3 induces no electric polarization, and so a straight- 
forward density-functional study of Cr 2 3 does not give 
information about the strength of the finite-temperature 
perpendicular magnetoelectric coupling. However, the 
tttl ordering shown in Fig. Q] (b) renders Cr 2 3 multi- 
ferroic and induces the electric polarization 

V = P z {mi) = 2 AS 2 . (4) 

We next extract A using ab initio methods by enforc- 
ing ttt-i- s P m ordering and calculating the magnetically- 
induced polarization V . 

First-principles calculations of the magnetoelectric cou- 
pling: We compute A and the spin exchange parame- 
ters using plane-wave density-functional theory, as im- 
plemented in the Vienna Ab-initio Simulation Pack- 
age (VASP)j2£ We use PAW potentials 21 for core- valence 
partitioning, and the local-spin-density approximation 
with a rotationally invariant Hubbard-C/ (LSDA+J7) for 
the exchange-correlation potential^ Our Hubbard U — 
2.0 eV is the same value that was taken for computing 
the magnetoelectric response at zero temperature i 1 ^ We 
deliberately do not include spin-orbit coupling, and all 
calculations are for collinear spin densities, so that A 
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corresponds only to polarizations induced by exchange 
striction. 

We find that a plane- wave cutoff of 500 eV and 
Monkhorst-Pack 2 ^ fc-point sampling of 4 x 4 x 4 are suffi- 
cient for computing the properties of interest. Note that 
much finer k meshes would be required for accurately 
resolving the magnetocrystalline easy axis^ but non- 
relativistic properties are well converged with our chosen 
parameters. 

We work with space group i?3c at the experimental 
volume 23 of 96. OA 3 and rhombohedral angle of 55.f3°. 
The internal coordinates are relaxed within our den- 
sity functional calculations for the tit4- magnetic con- 
figuration, yielding coordinates x — 0.f536 for Cr, and 
x = 0.9426 for O in Wyckoff positions 4c and 6e respec- 
tively. Subsequently, the Heisenberg exchange couplings 
Ji — J5,— corresponding to Cr-Cr distances of 2.65 to 
4.f A, are computed by fitting a Heisenberg Hamiltonian 
to DFT total energies of twelve different spin configura- 
tions with fixed ion coordinates in the hexagonal setting 
of Ric. This method is analogous to that employed by 

Finally, we compute A by enforcing the spin configu- 
ration of ttti an d re-relaxing the ionic coordinates in 
the rhombohedral unit cell. The resulting ionic config- 
uration has a polar lattice distortion. We compute the 
magnitude of V = 0.585 /iC/cm 2 using the Berry phase 
approach^ which allows us to extract A. We note that 
V is of the same order of magnitude as the polarization 
induced by exchange interactions in multifcrroics with 
collinear spins 

Monte Carlo simulations: Using Eq.Q we can express 
the magnetoelectric coefficient a|| in terms of spin corre- 
lation functions: 



d{P z ) 



dH z 



k B T 



(5) 
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where (. . .) denotes the thermal average at temperature 
T, ks is the Boltzmann constant and \xb is the Bohr 
magneton. 

In the mean-field approximation, equivalent to replac- 
ing the scalar products of spins Sq^-S^ by (S a j)-Sp.k + 
S QJ • (S^fc) in the expression for P z , one obtain o 16 ' 17 
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(6) 



where vq is the unit cell volume, in agreement with the 
simple argument given above. The mean-field expression 
qualitatively explains the observed temperature depen- 
dence of a || : It first grows, together with the order pa- 
rameter G z , as the temperature drops below Tjv, then 
subsequently decreases and vanishes at T = 0, together 
with the longitudinal magnetic susceptibility \\\- 

In Reffl8l an attempt was made to take into account 
the effects of spin fluctuations using a higher-order decou- 
pling scheme. This approximation fails, however, close to 
the transition temperature where spin fluctuations are 



large. We include spin fluctuations by calculating a» nu- 
merically using Monte Carlo simulations of a system of 
864 classical spins with exchange constants and the mag- 
netoelectric coupling A obtained from our first principles 
calculations, as described above. 

Since the correlation function in the right-hand side of 
Eq. ([5]) is zero unless a single antiferromagnetic domain is 
selected, we apply to our finite spin system a weak stag- 
gered field along the z axis, h(—) a , where a = 1,2,3,4 
labels magnetic sublattices. This mimics the easy mag- 
netic anisotropy of Cr2C>3 as well as the cooling in elec- 
tric and magnetic fields, used in measurements of mag- 
netoelectric coefficients to select the domain with a given 
sign of magnetic order parameter. The field strength, 
h = 0.165 meV, was chosen so that it is small compared 
to the scale of exchange interactions, but large enough to 
make the Monte Carlo results independent of h. 

Figure [2] shows the temperature dependence of a\\ ob- 
tained from Monte Carlo simulations (blue circles) and in 
the mean-field approximation described above (red solid 
line). The onset of the magnetoelectric response in our 
Monte Carlo simulations, as well as the sharp peak in 
the specific heat (inset), shows that the antiferromagnetic 
order sets in at ~ 290 K, close to the experimentally ob- 
served transition temperature Tjy = 307 K. The maximal 
value of the magnetoelectric coefficient obtained from our 
simulations is 0.9 x I0~ 4 (in CGS units), in excellent 
agreement with the experimental value of 1 .0 x 10~ 4 (see 
Refl28h. The maximum value is reached at ~ 240 K that 
compares well to T max ~ 260 K found in experiment. The 
mean field transition temperature (425 K) and maximal 
«H are significantly higher than the Monte Carlo values, 
indicating the importance of spin fluctuations in this ma- 
terial. 

Conclusions: We have presented the first ah initio cal- 
culations of temperature-dependent linear magnetoelec- 
tric responses. The quantitative agreement of our results 
with experimental data on Cr203 demonstrates that the 
dominant parallel magnetoelectric coupling in this ma- 
terial originates from non-relativistic exchange interac- 
tions between electrons. The strong temperature depen- 
dence of the magnetoelectric coefficient a|| underscores 
the general importance of spin fluctuations for magne- 
toelectric responses of materials with collinear spin or- 
ders. The magnetic-field-induced electric polarization 
of Cr203 is comparable to that predicted recently for 
an exchange-interaction-driven Kagome antiferromagnet 
with non-collinear spin ordering,! 1 However, the magne- 
toelectric response of non-collinear magnets does not van- 
ish at zero temperature as it involves transverse rather 
than longitudinal magnetic susceptibility and originates 
from the dynamic 'electromagnon' modes, which can be 
excited by both electric and magnetic fields and which 
are absent in collinear magnets i 29 ' 30 

The approach used in this paper, specifically the com- 
bination of first principles calculations for artificially im- 
posed magnetic states to extract parameters with Monte 
Carlo simulations of physically interesting quantities, 
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FIG. 2: (Color online) Temperature dependence of the mag- 
netoelectric coupling an obtained using ah initio values of the 
exchange constants and magnetoelectric coupling combined 
with Monte Carlo simulations (blue circles) and mean field 
calculations (solid red line). The inset shows the temperature 
dependence of magnetic specific heat. 



opens a route to theoretical studies of a large variety of 
temperature-dependent static and dynamic magnetoelec- 
tric phenomena. Accurate predictions of the magnitude 
of magnetoelectric responses at finite temperature will 
greatly facilitate the search for and design of materials 
with the strongest responses. 
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